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Condensed  Storage  of  Diffusion  Equation  Solutions  for  Atmospheric  Density 
Model  Computations. 

1.  INTRODUCTION 

1  2 

Due  to  the  increasing  complexity  of  recent  Jacchia  density  models  *  it 

becomes  necessary  to  store  in  tabular  form  the  solution  of  the  diffusion 

equation  for  each  constituent,  even  if  only  the  total  mass  density  is  to 

be  computed.  This  is  because  the  models  call  for  corrections  to  the 

diffusion  process  for  such  effects  as  local  time  phase  variations, 

seasonal-latitudinal  variations  and  geomagnetic  activity,  all  of  which 

differ  for  the  different  constituents.  The  Jacchia  1964  model"1,  on 

the  other  hand,  required  only  temperature  corrections,  and  hence  only  the 

total  mass  density  needed  to  be  stored  as  a  function  of  height  and 

exospheric  temperature.  As  a  result,  one  must  either  set  aside  a  large 

2 

block  of  storage  (the  current  SAO  version  of  the  Jacchia  1977  model 
would  require  in  excess  of  17,000  words),  or  resort  to  random  access  disk 
storage.  The  first  approach  is  undesirable  in  almost  all  cases, 
particularly  with  small  computers.  The  second  poses  problems  in 
implementation  on  different  computers. 

A  third  approach  is  presented  here  which  takes  advantage  of  the  fact  that, 
with  only  two  exceptions,  all  terms  in  the  diffusion  equation  are 
analytically  integrable.  The  exceptions  are  the  scale  height  term,  which 
differs  for  different  constituents  merely  by  a  mass  factor,  and  the 
vertical  flux  term  for  hydrogen,  which  is  important  only  for  low 
altitudes.  It  therefore  is  necessary  to  store  only  the  integral  of  the 
corrmon  mass- independent  factor  of  the  former  term,  and  the  solution  for 
hydrogen  whenever  that  constituent  is  of  special  interest. 
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2. 


ANALYSIS 


The  diffusion  equation  for  the  ith  constituent,  in  the 
(J77)  density  model  is2*** 


dn^  dT 

n .  +  T 

1 

where 


(1 


+ai5 


dh 
n . 


0 


=  ith  constituent  number  density 
T  =  temperature 

'» ^  =  ith  constituent  thermal  diffusion  coefficient 
h  =  height 
H1  =  R*T  /M-j  g 

R  =  8.31432X103  J(kg  -  mol)"1/* 

=  ith  constituent  molecular  mass 
g  =  9.80665  (1  +h/Re)"2m/sec2 
Re  =  6.356766X106m 

=  ith  constituent  vertical  flux 
D  =  diffusion  coefficient  =  2X1020  f T/N 
N  =  total  number  density 


The  are  assumed  values  of  -0.38  and  -0.25  for  helium 
and  0  for  all  others.  The  41  are  0  for  all  but  hydrogen. 

*  j 

the  vertical  flux  term  leads  to  * 


«i  U 

where 


"i  <V  u 


K' 

exn  f 

TTh,  TJ 

9  L 

[H1  F(h, 


F(h  ,  h,  TJ 


fh  h  g(z)/R*T{z ,  TJ  d z 


T^=  exospheric  temperature 
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Jacchia  1977 


and  hydrogen 
Neglecting 


2.1 


Tabulation 

In  practice  it  is  therefore  necessary  to  tabulate  only  F  and  the  density 
for  one  of  the  constituents  at  hQ  as  a  function  of  Tt  ,  if  hQ  is  chosen 
to  be  the  homopause,  100km.  If  ^  is  chosen,  then  the  others,  except 
hydrogen,  are  given  by 

log  ni  (hQ,  T  ,  )  =  log  n^  ( hQ ,  T  )  ♦  Qi 


where  the  subscript  r.  ■  .ates  species  and  the  are  constants: 

q,  -  log  (q,/q28)  ’  '  lc'  32 

Qie  =  -log  q2s  -log  (M,/R0‘) 

+  log  L2  (1-PC /'Nq1  )] 


q32  =  -log  g28 


-log 
+  lon 


(M '/*</) 


(i  +  q 


Where 


=  sea  level  concentration  of  i 


th 


constituent 


M*  =  mean  molecular  mass  at  100km  (Eq.  5,  ref.  2) 

M  1  =  mean  molecular  mass  at  sea  level 

o 

If  one  is  interested  in  including  the  escape  flux  term  for  hydrogen, 
special  tables  would  still  be  necessary;  however  the  storage  requirement 
would  still  be  considerably  less  than  if  separate  tables  are  used  for  all 
constituents.  Furthermore,  since  the  escape  flux  term  is  important  only 
below  500km,  it  would  be  necessary  to  store  results  only  for  that  region. 
If  one  excludes  the  escape  flux  term,  the  H  density  is  computed  using: 
hQ  =  500km;  eq.  17  of  ref.  2: 
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log  nj  (500,  )  =  5.94  +  28.9  1/4, 

and 

F  (500,  h,  T:,  )  «  F  (100,  h,  T„  )  -  F  (100,  500,  T  ). 

2.1.1  Homogeneous  Layer  (90km  <  h  <  100km) 

For  the  homogeneous  layer  the  diffusion  equations  are  replaced  by  a 
single  barometric  equation  for  the  mass  density,  from  which  component 
densities  may  be  derived,  as  indicated  in  reference  2.  Hence  only  the 
mass  density  need  be  stored.  Alternatively  one  may  store  the  density  for 
and  derive  the  others  from  it,  as  for  the  homopause  boundary.  The 
equations  are  the  same  except  that  M‘  would  be  the  mean  molecular 
weight  at  the  height  of  interest,  given  by  Eq.  5.  of  ref.  2. 


4 


l 


3. 


RESULTS 


A  subroutine  has  been  constructed  to  implement  the  above  procedures. 

Atomic  nitrogen  has  been  included,  based  on  the  most  recent  SAO 

subroutine  version,  which  incorporates  the  AE  OSS  mass  spectrometer 
5  6 

data.  ’  For  this  we  have  taken  the  fractional  sea-level  composition 
given  by  Table  1,  with  the  resulting  mean-mol  ecu!  ar  weight  Mo  = 

'  =  28.9586. 

The  tables  for  the  diffusion  function  F  are  built-in  along  with  the  N~ 
densities  for  the  homogeneous  layer  and  the  atomic  hydrogen  solution 
including  escape  flux.  Together  these  occupy  less  storage  than  required 
for  a  single  height  segment  (4000  words)  when  the  solutions  are  tabulated 
for  each  component  on  a  random  access  file.  Hence  no  external  files  are 
required.  Run  times  on  the  CDC  6600  are  comparable  to  the  random  access 
version,  with  essentially  identical  results. 

A  FORTRAN  deck  of  this  subroutine  is  available  on  request. 


Table  1. 


Sea-Level  Composition 


Constituent 


Fraction  by  Volume  Molecular  Mass 


Molecular  Nitrogen  (N2) 
Oxygen  (0^) 

Argon  (Ar) 

Hel lum  (He) 

Atomic  Nitrogen  (N) 


0.78103 

28.0134 

0.20553 

31.9988 

0.009342 

39.948 

0.000005242 

4.0026 

0.00007502 

14.0067 
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